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ABSTRACT 


Zero-field  parabolic  states  of  hydrogen  are  used  as  a  basis  set  for  the  study  of  low-order 
harmonic  generation  in  moderately  strong  laser  fields.  The  disappearance  of  bound  parabolic 
states  with  large  electric  dipole  moments  in  moderately  strong  fields  leads  to  the  simplification  of 
the  expression  for  the  total  time-dependent  dipole  moment  of  the  atom.  By  including  only 
contributions  from  bound  atomic  states  and  with  some  simplifying  assumptions  it  is  possible  to 
express  the  total  dipole  moment  by  an  analytic  formula.  The  present  results  are  compared  to  the 
benchmark  calculations  [1,2]. 
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1.  INTRODUCTION 

The  interaction  of  an  intense  ultrashort  pulse  laser  with  an  atom,  molecule,  cluster,  or 
solid  target  gives  rise  to  a  number  of  competing  processes  including  the  generation  of  harmonic 
radiation  [3].  The  major  focus  on  harmonic  generation  has  been  on  the  successful  generation  of 
higher  order  harmonics  due  to  the  interest  in  developing  shorter  wavelength  tunable  coherent 
radiation  sources  in  the  VUV-  and  X-ray  regimes  of  the  spectrum  for  applications  in,  for 
example,  plasma  diagnostics,  lithography,  and  x-ray  holography.  The  spectral  distribution  of  the 
harmonic  intensity  as  a  function  of  harmonic  number  exhibits  a  characteristic  spectra  displaying  a 
rapid  decay  of  intensity  for  low  harmonics  followed  by  a  plateau  region  where  the  intensity  is 
relatively  insensitive  to  harmonic  number  and  remains  fairly  flat  only  finally  to  fall  precipitously 
for  higher  harmonics.  The  shape  of  the  spectra  for  the  plateau  and  cutoff  harmonics  can  be 
reasonably  explained  in  terms  of  a  semiclassical  recollision  model  often  referred  to  as  the 
tunneling- recombination  model  that  predicts  a  cutoff  energy  of  -f  3.1717^  where  is  the 

ionization  energy  and  U ^  is  the  ponderomotive  energy  of  the  laser  field.  Once  the  electron  is 

liberated  from  the  core  it  is  accelerated  in  the  field  until  the  phase  of  the  field  reverses  and 
pushes  the  electron  back  toward  an  ion  where  it  may  radiatively  recombine  leading  to  harmonic 
generation  [4,5].  The  semiclassical  model  has  been  useful  for  explaining  some  of  the 
experimentally  observed  features.  Since  its  inception  the  semiclassical  model  has  become  the 
standard  with  which  a  number  of  classical,  quantum,  and  quantum  field  theoretic  models  have 
been  compared.  It  is  not  our  intent  to  discuss  the  pros  and  cons  of  each  model  and  therefore 
refer  the  interested  reader  to  three  excellent  review  articles  [6-9]  where  they  compare,  contrast, 
and  consider  the  advantages  and  disadvantages  of  each  approach.  The  most  recent  quantum 
field  model  appears  in  [10- 13]  and  contain  adequate  references  to  earlier  works.  However,  by 
their  very  nature  many  of  these  models  are  very  complicated  and  computer  intensive  and 
although  they  provide  insight  and  explanations  for  the  purely  quantum  effects  they  lose  the 
simplicity  characteristic  of  the  semiclassical  model.  Therefore,  in  the  spirit  of  the  semiclassical 
model  we  describe  what  we  believe  to  be  a  simple  theory  that  includes  some  of  the  effects  of 
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the  quantum  and  classical  models  but  remains  physically  transparent  and  calculationally 
straightforward. 

The  physics  responsible  for  harmonic  generation  is  due  to  the  asymmetry  in  the 
distribution  of  the  electronic  charge  density  in  each  atom  composing  the  medium  by  the  time- 
varying  driving  laser  field.  The  action  of  the  field  produces  a  time- dependent  dipole  moment 
whose  phase  is  related  to  the  incident  applied  laser  field  and  whose  Fourier  transform  is  related 
to  the  strengths  and  phases  of  the  harmonic  radiation  from  the  individual  atoms. 

It  was  not  known  initially  if  the  harmonics  were  dominated  by  single  atom  interactions  or 
if  propagation  through  the  medium  controlled  the  spectral  behavior.  Actually,  single  atom 
calculations,  even  simple  two- level  atom  approaches,  have  reproduced  many  of  the  features 
associated  with  strongly  driven  nonlinear  systems.  The  characteristic  profile  of  harmonic  order 
as  a  fiinction  of  dipole  moment,  including  the  plateau  and  cutoff,  is  often  explained  in  terms  of 
field  induced  ionization  that  becomes  the  source  term  for  electrons  bom  with  zero  initial  velocity 
in  the  laser  focus.  Following  liberation  fi'om  the  atom  some  of  the  electrons  follow  trajectories 
on  the  laser  half-cycle  that  return  to  recollide  and  radiatively  recombine  with  the  nucleus 
producing  harmonic  radiation.  This  semiclassical  picture  captures  many  of  the  observed  features 
of  harmonic  generation  including  the  correct  scaling  of  the  harmonic  spectrum  with  the  intensity 
of  the  field.  In  general,  the  modeling  of  strong- field  harmonic  conversion  requires  both  the 
evaluation  of  the  single  atom  response  and  the  solution  of  the  propagation  equations  for  the 
harmonic  fields.  The  propagation  calculation  requires  knowledge  of  the  single-atom  dipole 
strength  (and  phase)  as  a  function  of  the  intensity  of  the  driving  field,  for  a  large  number  of 
intensity  points,  typically  on  the  order  of  hundreds.  This  explains  why  most  calculations  have 
been  carried  out  using  simple  models  for  describing  harmonic  generation.  In  this  investigation  we 
will  focus  on  a  simplified  approach  for  the  single  atom  aspects  of  the  problem  and  refer  to  the 
work  of  Borisov  et  al.  [14]  for  the  propagation  issues. 

It  is  also  worth  mentioning  that  the  efficiency  of  high  harmonic  generation  is  not  yet  high 
enough  for  many  of  the  potential  applications  mentioned  above.  The  problem  is  twofold;  1)  not 
enough  energy  can  be  stored  in  the  higher  harmonics  and  2)  phase  mismatch  between  the 
fundamental  and  the  harmonic  fields  can  result  in  the  amplitude  of  a  given  harmonic  changing 
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dramatically  with  phase  upwards  of  3  orders  of  magnitude  depending  on  the  specific  case.  A 
possible  resolution  of  this  problem  might  be  accomplished  by  channeling  the  emitted  radiation 
into  a  single  harmonic  by  the  process  of  pulse  shaping  [15].  Finally,  since  much  of  the  recent 
emphasis  has  been  on  higher-order  harmonics  we  focus  on  the  low-order  harmonics  and 
demonstrate  that  they  can  be  described  adequately  by  the  bound  atomic  states. 

II.  MODEL  AND  CALCULATIONS 

There  are  basically  two  quantum  mechanical  approaches  to  the  theoretical  study  of 
harmonic  generation;  the  first  one  is  based  on  the  direct  numerical  solution  of  the  time- 
dependent  Schrodinger  equation  [1,16-18],  the  other  one  on  the  Floquet  formalism  [2,19,20]. 
If  accurate  results  are  required,  both  approaches  are  computationally  demanding.  This  has  led 
to  a  number  of  modified  and  simplified  treatments  of  harmonic  generation  with  various  degrees 
of  complexity.  Some  procedures,  besides  being  simple,  may  also  offer  certain  insights  into 
individual  physical  processes  that  contribute  to  the  generation  of  harmonic  radiation. 

The  recollision  model  of  harmonic  generation,  even  in  its  fully  quantum  theoretical  form 
[21,22],  does  not  include  aiy  participation  of  excited  bound  states  in  the  emission  of  high 
harmonics  and  is  essentially  valid  for  harmonics  with  energy  well  above  the  ionization  energy 
[21].  On  the  other  hand,  several  previous  studies  [23-29]  showed  that  excited  levels  may  be  of 
inportance  for  harmonic  generation.  In  the  present  paper  we  investigate  the  effect  of  radiation 
on  excited  levels  and  show  that  the  harmonic  emission  of  hydrogen  at  moderate  laser  intensities 
may  be  explained  by  transitions  between  the  side-bands  of  these  levels  and  the  ground  level. 

We  study  the  hydrogen  atom  in  the  presence  of  a  moderately  intense  linearly  polarized 
radiation  field  of  infinite  duration  and  look  for  an  approximate  stationary  solution  of  the 
Schrodinger  equation  in  order  to  obtain  a  simple  analytic  expression  describing  the  harmonic 
radiation. 

The  wave  function  of  a  hydrogen  atom  in  a  linearly  polarized  radiation  field  may  be 
conveniently  expressed  in  terms  of  zero- field  parabolic  states  |«  n,  w)  .  An  important  property 
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of  these  states  is  the  generally  asymmetric  charge  distribution  with  respect  to  the  xy  plane, 
indicating  an  increasing  possibility  of  field  ionization  for  some  states  with  nonvanishing  values  of 
the  dipole  moment  (n  «,  m  \z\n  «,  n^.  In  order  to  obtain  a  simple  expression  for  the  dipole 

moment  of  the  atom  in  the  radiation  field,  including  only  contributions  fi'om  bound  states,  we 
adopt  the  following  two  simplifications:  the  continuum  states  of  the  atom  will  be  ignored  and  all 
matrix  elements  {n  «,  m  \z\n'  n\  m'^  taken  between  excited  states  with  n^n'  will  be  omitted. 

These  simplifications  allow  us  to  express  the  total  induced  dipole  moment  of  the  atom  in  closed 
form  as  a  function  of  the  field  amplitude  and  a  set  of  atomic  parameters.  With  the  electric  field 
in  the  z- direction,  the  part  of  the  Hamiltonian  representing  the  interaction  between  the  bound 
electron  and  the  field  in  the  dipole  approximation  has  the  form  Fzsin®  / .  In  our  procedure 
the  \n  m)  states  are  coupled  only  to  the  ground  state  and  the  total  wave  function  'F  will 

contain  only  parabolic  states  with  m  =  0.  In  the  subsequent  equations  the  label  m  will  be 
omitted  and  atomic  units  wilt  be  used  throughout  the  paper  unless  indicated  otherwise. 

In  the  basis  of  parabolic  states,  the  matrix  with  elements  (n  n^  \z\  n  n,)  is  diagonal  and 


(nn^  I  z  I  ««|)  =  ~  ”2 ) 


(1) 


with  «2  =n-n^-\.  Factors  {n^-n^)  assume  values  -(«-l), -(n-3),...«-3,  n-1 .  For 
our  purpose,  it  is  convenient  to  designate  each  parabolic  state  by  |/7  A )  according  to  the  value 
of  the  matrix  element  (l)  so  that 
^nx  ~  I  z  I  n  A)  =  3n  A  for  n  odd, 


=  {nX  I  z  I  nX)  =  3n( A  ±  — )  for  n  even  and  A  0. 

2  > 


(2) 


For  n  odd,  A  =0,±1, ±2,. .+—(«-!),  and  for  n  even,  A  =  ±1, ±2,- •  ;±—« .  The  functions 

2  2 

|n  A  )  of  the  zero- field  parabolic  states  are  superpositions  of  regular  spherical  functions  (p^ 
represented  as 

(3) 
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and  can  be  expressed  in  terms  of  S-j  symbols  [30,31].  With  the  convention  that  are 
positive  for  f = 0 ,  we  have 


^(n-l)/2  (n-t)/2 


a 


-a  0 


(4) 


where  a  =  X  for  «  odd,  X  ±  for  «  even  and  X  0 ,  with  special  values 


> 


Following  the  Floquet  method,  we  require  the  total  wave  fimction  to  have  the  form 


=X  I”^)  exp[-/(£’,  +  e)/]. 

nX 


(6) 


is  the  energy  of  the  zero- field  ground  state  and  (£'i  +  e)  is  the  quasi-energy  of  the 


system.  Due  to  the  periodicity  of  the  Hamiltonian  it  is  possible  to  write  {t)  in  the  form 

Kx  (0  =  C„x  (Oexp (/ cos  cot)  (7) 

and  the  solution  of  the  Schrddinger  equation  may  be  obtained  in  a  convenient  form  if  constants 
are  properly  chosen.  With  the  definition 

(8) 

the  matrix  elements  of  the  operator 

i- — {H^  +  Fz  sinmt), 
at 


where  is  the  field- free  Hamiltonian,  are  zero  in  the  basis  of  fimctions 

I nX) exp(/ cos  cot—iE^ t) 

for  a  given  n .  Therefore  the  strongest  effect  of  the  field  on  the  atom  is  already  properly 
included  in  the  basis  set  if5„;^  in  (6)  are  defined  according  to  (7)  and  (8).  Parabolic  states  |«0) 

are  not  coupled  to  the  ground  state  |l0)  due  to  (5)  and  they  wiU  not  be  included  in  the 
expansion  (6).  Using  the  substitution 

C.o  =  =  (9) 

p  ? 
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together  with  the  expansion  of  the  exponential  factor  in  (7)  in  terms  of  Bessel  functions 
J  derives  a  system  of  coupled  equations  for  and  (with  «  >  1).  We 

require  a  solution  for  which  =  1  in  the  limit  of  vanishing  field.  Then  fi'om  the  properties  of 
and  it  follows  that  all  C,(,  Avith  p  odd  should  be  zero  and  that  the  matrix  of 
coefficients  for  and  is  real  and  symmetric.  Moreover  we  find  that 

c;?=(-irc“.  (10) 

The  sum  over  two  terms  with  opposite  values  of  A  may  be  readily  evaluated  and  the  system 
may  be  further  reduced  to  the  positive  values  of  A  . 

With  A  >  0  and  p  even,  we  obtain 

(e+£,-«.-?®)4”+5:^j„c“  =  o,  .  (11) 

P 

Kxp,  =  Fa[\0,n\)Ku{p-^)tlJp-,[Lx)- 

a{ni,n' a  matrix  element  <  >  taken  between  spherical  states  with  w=0,  and 

are  zero  field  energies  of  the  kA  states.  The  summation  over  p  and  q  is  fi'om  -o®  to  +oo  and 
the  summation  over  nX  includes  all  parabolic  states  \nX)  for  n  >  1  and  A  >  0  which  may  exist 

at  a  particular  amplitude  F  of  the  electric  field.  Conditions  for  the  inclusion  of  \nX)  states  to  the 

sum  in  (11)  will  be  discussed  later  in  par.  in.  The  solution  of  (11)  yields  relative  values  of 
coefficients  and  eigenvalues  e. 

To  obtain  the  dipole  moment  of  the  atom  in  the  radiation  field,  we  chose  of  all  possible 
solutions  of  (1 1)  the  one  with  the  largest  value  of  ,  which  has  a  character  of  the  unperturbed 
ground  state.  Such  a  solution  corresponds  to  the  smallest  absolute  value  of  e,  at  least  for  small 
values  of  the  field  amplitude  F .  The  time- dependent  dipole  moment  in  the  length  gauge  is  given 
by 
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(12) 


Z)(/)  =  {'P|z|4') 

and  it  may  be  written  as  a  sum  over  terms  with  time  factors  ,  in  evaluating 

(12)  we  include  only  terms  vrith  matrix  elements  connecting  the  ground  level  with  the  excited 
nX  levels  which  are  expected  to  give  dominant  contributions  to  the  dipole  moment.  Matrix 
elements  ^10|z|/7A^are  real  and  the  summation  over  ±X  in  the  expression  for  E>{t)  leads  to 

cancellation  of  all  terms  with  factors  for  which  q +11- p  is  even.  We  define  an  odd 

integer 

N  =  \q  +  p-p\  (13) 

and  obtain 

D(?)  =  E(£)^e'''"'  +  c.c.)  (14) 

with 

D,  =2  I  05) 

nXpq^ 

As  in  (1 1),  the  summation  over  parabolic  states  \nX)  in  (15)  is  limited  due  to  the  presence  of 
the  electric  field,  p  is  even,  A  >  0  and  p,  q,  p  may  have  values  fi'om  — to  +o°  and  should 
satisfy  the  condition  (13).  Amplitudes  C|^^and  in  (15)  should  be  normalized.  The  rate  of 
emission  of  harmonic  photons  with  energy  (o'  =  Nco  by  an  isolated  atom  is  then 

(16) 

III.  EFFECTS  OF  IONIZATION 

The  simplicity  of  expression  (15)  for  is  achieved  mainly  because  of  the  omission  of 
continuum  states  in  expansion  (6)  for  the  total  wave  function  and  to  factorization  of  in  (7). 
An  important  consequence  of  our  simplifying  assumption  is  the  neglect  of  all  processes  leading 
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to  ioni2ation.  Field  ionization  of  excited  nX  levels  will  decrease  the  availability  of  these  levels  for 
the  generation  of  harmonic  emission.  Maintaining  the  simplicity  of  the  expression  forD^^, ,  one 
may  approximately  take  into  account  the  effects  of  ionization  due  to  the  electric  field  by 
multiplying  all  contributions  fiom  nX  levels  to  the  dipole  moment  (15)  by  a  factor  representing 
the  expectation  value  of  these  levels  in  the  AC  field.  To  determine  for  a  given 

amplitude  of  the  field,  we  will  assume  first  of  all  that  P^^  is  zero  if  F  is  substantially  larger  than 
a  certain  critical  amplitude  .  In  the  static  field,  the  threshold  value  of  the  field  strength 
necessary  for  ionization  of  the  n  -level  is  (2n)‘^ .  However,  for  parabolic  states  F^  will  depend 

on  X  within  the  same  n  -manifold  due  to  the  asymmetry  of  the  charge  distribution  with  respect  to 
the  xy-plane.  This  effect  has  been  studied  for  static  fields  by  Banks  and  Leopold  [32],  Harmin 
[33],  and  Popov,  et.  al.  [34],  A  better  approximation  applicable  to  parabolic  states  with  largest 
tunneling  ionization  rate  for  a  given  n  and  field  in  the  +z  direction,  i.e.  with 

A  =  — ^—n  for  n  even  and 

2 

A  =  — -  1)  for  n  odd,  yields 


F^  =0.130n  ^  (Krainov,  et.  al.  [35])  (17) 

For  such  field  strength  the  distance  of  the  saddle  point  of  the  potential  fi-om  the  nucleus  along  the 
z  -axis  is  equal  to  r,  =F~^''^  ~  l.lln^ .  This  expression  is  similar  to  the  expression  for  \d„j\  of 
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these  levels,  namely  |t/„;i|  =  — for  «  even  and  — «(«-!)  for  n  odd.  The  similarity  becomes 


even  more  apparent  if  we  consider  that  (17)  represents  only  the  lower  limit  of  the  correct  value 
of  F^  [36].  We  may  therefore  approximately  find  the  value  of  the  critical  field  by  setting 

Fc-d:l.  (18) 


According  to  eq.  (18)  the  center  of  gravity  of  the  electron  mass  coincides  with  the  saddle  point 
of  the  potential.  This  is  in  agreement  with  the  fact  that  ionization  can  really  occur  just  in  or  near 
the  region  where  the  Coulomb  field  is  smaller  than  the  external  field,  and  therefore  we  may 
generalize  the  validity  of  (1 8)  to  all  A  <  0  states  of  the  same  n  -manifold  and  to  field  strengths 
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even  larger  than  (In)"'* .  Eq.  (18)  is  consistent  with  the  strong  X  -  dependence  of  ionization 
rates  found  by  Luc-Koenig  and  Bachelier  [37]  in  a  static  electric  field. 

In  the  AC  field,  the  A  -  dependence  of  is  due  to  the  splitting  of  the  n  -levels  into  the 

n£  or  equivalently  to  nX  components  (Delone  and  Krainov  [38])  and  to  overlap  of  these 
components  with  adjacent  «  +  !  levels.  Assuming  the  validity  of  (18)  also  for  the  oscillating 
fields  and  for  levels  with  A  >  0 ,  then  for  the  most  weakly  bound  parabolic  states  with  a  given 
n ,  behaves  as  ,  and  for  the  most  tightly  bound  states  (with  smallest  non- vanishing 

Kl).  F^  scales  as  «  ^  in  the  limit  of  large  n . 

Our  calculations  were  performed  for  radiation  intensities  fi'om  lO'^  to  8x1 0'^  W/cm^  for 
the  wavelength  1064  nm  (to  =  0.04282).  The  lowest  intensity  corresponds  to  F  =0.005345 
and  to  1=13.68  according  to  (18).  At  this  intensity  only  the  following  excited  nA  levels 
satisfy  the  condition  \  |<13.68; 

n  |A|  I^J 

2  1  3 

3  1  9 

4  1  6 

6  1  9 

8  1  12 

We  note  that  only  levels  with  the  smallest  value  of  |  for  each  n  are  present  above 
F  =0.005345  and  thus  the  summation  over  nX  in  (15)  reduces  to  a  very  small  number  of 
terms. 

At  the  intensity  of  8xl0'^  W/cm^  only  the  lowest  excited  state  with  «  =  2  could  be 
present  according  to  the  condition  (18).  In  reality  the  expectation  value  Kx  continuously 

changes  in  a  broad  range  around  F^  and  it  depends  on  the  actual  process  of  ionization  and 
therefore  on  the  Keldysh  parameter  y.  This  parameter  varies  fi'om  4.0  to  0.1 1  fiom  n  =  2  to  n 
=  8  in  the  entire  interval  of  intensities.  Consequently,  no  unique  process  will  be  solely 
responsible  for  destruction  of  bound  states  at  all  intensities.  For  simplicity  we  have  adopted  the 
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photoionization  as  the  dominant  process  and  determined  PnX  for  all  considered  intensities  from 
the  expression 

P„,  =  {\  +  2nn^W^r  (19) 

where  2%  « ^  is  the  Kepler  orbital  period  for  the  « th  level  and  is  the  multiphoton  ionization 
rate  of  the  nX  state.  Because  of  the  approximate  character  of  the  expression  (18)  for  the 
critical  field  ,  we  have  also  included  levels  with  n  =5  (|  |=1 5). 

Multiphoton  ionization  rates  W  were  calculated  using  a  Keldysh- like  method  in  which 
the  final  state  is  represented  by  a  Coulomb-Volkov  continuum.  This  method  was  proposed  by 
Pan  [39]  and  later  used  in  ref  [40],  where  a  detailed  description  is  given.  Because  of  the  value 
of  (D,  levels  with  n  =2  and  3  need  more  than  one  photon  to  be  ionized.  Table  1  shows  values  of 
for  all  included  levels  and  intensities  of  radiation. 

Table  1. 

Expectation  values  P„x  of  parabolic  states  n,\X\  =  \. 


2 

3 

4 

5 

6 

8 

1(1  o'\ 

W/cm^)\ 

0.1 

0.972 

0.958 

0.620 

0.671 

0.761 

0.728 

0.5 

0.628 

0.482 

0.243 

0.322 

0.391 

0.359 

1 

0.411 

0.212 

0.138 

0.180 

0.248 

0.233 

2 

0.357 

0.112 

0.079 

0.125 

0.158 

0.161 

4 

0.314 

0.157 

0.054 

0.115 

0.121 

0.133 

8 

0.274 

0.206 

0.044 

0.105 

0.109 

0.129 

W  e  note  that  P^^^  decreases  more  rapidly  for  n  =2,3  than  for  other  levels,  obviously  as  a 
consequence  of  non-existent  single  photon  ionization  for  these  levels.  The  decrease  of  P^j^  for 
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all  levels  becomes  slower  at  intensities  above  lO'^  W/cm^,  and  the  «=3  level  shows  an 
indication  of  a  possible  onset  of  stabilization. 

IV.  RESULTS  AND  DISCUSSION 

The  results  of  our  calculation  for  hydrogen  according  to  (15)  with  the  inclusion  of 
factors  from  Table  1  are  shown  on  Figs.  1  and  2  for  wavelength  =  1064  nm  (to  =  0.04282) 

and  I  =  10‘^  to  8xl0'^  W/cm^.  The  harmonic  spectra  on  the  Figures  are  compared  to  grid 
calculations  of  Krause  et.  al.  [1]  for  I  >  lO'^  W/cm^  and  to  results  of  the  Floquet  method  by 
Potvliege  and  Shakeshaft  [2]  for  1  <  lO'^  Relative  data  for  emission  rates  in  [2]  were  first 

converted  to  according  to  (16)  and  then  to  absolute  values  using  values  of 

calculated  by  Krause  et.  al.  [1].  This  appears  well  justified  due  to  excellent  agreement  between 
relative  values  of  both  results.  Data  ior  N  =1  were  taken  from  [2]  for  all  values  of  I. 

Under  certain  conditions  expression  (15)  for  may  be  further  simplified.  For  low 

values  of  F  the  coefficients  for  p  ^0  are  very  small  and  the  eigenvalue  e  approaches 

zero.  In  our  calculation  for  intensities  up  to  8x10'^  W/cm^  for  hydrogen  we  observed  that  the 
contribution  of  terms  with  p  were  small  and  that  |e|  was  always  smaller  than  0.05. 

Moreover,  if  Oi  is  such  that  [i?,  -  -  qoi\  is  not  close  to  zero  for  any  n  and  ^ ,  it  is  possible 

to  find  an  approximate  expression  for  from  the  second  equation  in  (11)  by  setting 
Cj''^  =  S(p,0),£  =  0 .  Substitution  into  (15)  and  inclusion  of  factors  leads  to  a  closed  form 
for  , 

D,  =  2i^Fj^[a{\0,n\)K,,f{E, -E„  -geo)-'  +J-N-AL,mx-  (20) 

nXq 

The  summation  over  nA  and  q  is  the  same  as  in  (11).  The  harmonic  spectrum  of  hydrogen 
according  to  (20)  agrees  in  most  cases  within  a  factor  of  2  with  results  presented  on  Figs.  1  and  2. 
The  difference  decreases  with  increasing  N  and  with  decreasing  intensity  of  radiation. 

The  basic  features  of  harmonic  spectra  are  well  reproduced  in  the  present  calculations: 

steep  decline  of  \Dfjf  for  low  harmonics,  the  onset  of  a  plateau  in  the  vicinity  of  N=1  and  then  a 
rapid  decrease.  This  behavior  and  the  value  of  the  plateau  may  be  explained  by  the  properties  of 
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Bessel  functions  as  discussed  in  detail  by  Burlon  et  al.  [27]  in  the  case  of  a  degenerate  two-level 
model  atom. 

It  has  been  shown  by  Krause  et.  al.  [1]  that  \D^  f  of  higher  harmonics  exhibits  rapid 
fluctuations  of  up  to  one  order  of  magnitude  as  a  function  of  the  radiation  intensity  I.  We  have 
found  similar  behavior  in  our  calculations,  but  various  local  maxima  and  minima  do  not  generally 
agree  with  [1]  so  that  occasionally  large  discrepancies  may  occur  for  a  particular  value  of  I  or  N  , 
e.g.  for  1  =  4xl0'^  W/cm^,  N  =11  and  I  =  8xl0’^  W/cm^,  N  =  3,19.  Nevertheless,  some  details 
of  the  I- dependence  of  several  harmonics  reported  in  [1]  are  well  reproduced  in  our  results,  e.g. 
the  broad  and  flat  maximum  between  I  =  10'^  and  2x10’^  W/cm^  of  =  7  and  9,  and  a  smooth 
behavior  of  =11  and  19  in  approximately  the  same  region. 

The  time- dependent  dipole  moment  of  the  atom  in  the  velocity  gauge  is  defined,  instead  of 
(12),  by 


„  V  iNo)  t 
=  L  +  cc). 

N 


(21) 


Because  (  10 


dz 


nX)  =  (E^-E^)a(\0,nl)b^^;^,D 


multiplying  each  term  in  (15)  by  a  factor  -/(£„  -E^). 


may  be  obtained  from  by 


The  emission  rate  F  of  harmonic  photons  with  energy  (O  =  N(0  is  then  again  given  by 


V  —1 

(16),  but  should  be  replaced  by  D^{Nco)  .  As  a  result,  the  emission  rate  obtained  from  the 
velocity  expression  is  close  to  the  rate  from  the  length  formula  for  such  harmonics  for  which  the 
ratio  (E^  ~E^)I  NCO  of  contributing  nX  levels  is  close  to  1,  i.e.  for  values  of  N  around  1 1. 

For  lower  N  ,  the  velocity  rate  is  higher,  and  for  larger  N  it  is  lower  than  the  length  rate.  For 
example,  at  8xl0'^  W/cm^  the  velocity  rate  is  approximately  10  times  larger  than  the  length  rate  for 
N  =3,  and  10  times  smaller  for  N  =33. 

For  the  intensity  of  8x10'^  W/cm^  and  even  greater  intensities,  higher  harmonics  require  the 
presence  of  levels  with  \  >  20,  but  such  levels  cannot  exist  according  to  our  assumption 
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regarding  the  critical  value  .  Therefore,  the  mechanism  of  harmonic  generation  in  hydrogen 
described  in  this  paper  is  questionable  beyond  10’''  W/cm^  unless  a  stabilization  of  high  bound 
levels  takes  place. 

In  conclusion,  it  was  shown  that  with  some  simplifying  assumptions  and  zero- field  bound 
parabolic  states  as  basis  set  it  is  possible  to  express  the  time- dependent  dipole  moment  of  a 
hydrogen  atom  in  a  radiation  field  by  a  simple  analytic  formula  (20).  According  to  the  present 
description,  the  harmonic  radiation  of  photons  with  energy  up  to  two  times  the  ionization  energy 
may  be  interpreted  as  the  result  of  transitions  between  the  ground  state  and  the  energy  side¬ 
bands  of  excited  states,  which  are  coupled  to  the  ground  level  by  the  incident  radiation.  The 
number  of  parabolic  states  included  in  the  expression  for  the  dipole  moment  depends  on  the 
radiation  intensity,  and  the  critical  amplitude  excluding  levels  with  higher  values  of  \d„^  \  is 
given  by  (18).  Consequently,  in  the  interval  of  radiation  intensities  beyond  lO'^  W  cm'^,  only 
the  parabolic  levels  with  |A|  =  1  may  exist  and  scales  as  .  According  to  the  present 

results,  almost  all  harmonics  in  the  investigated  intensity  range  may  be  explained  just  by  the 
contributions  from  bound  states.  Our  assumptions  restrict  the  applicability  of  this  procedure  to 
radiation  intensities  less  than  10'“^  W/cm?  for  hydrogen. 
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CAPTIONS  TO  FIGURES 


Figure  1.  Harmonic  intensities  of  hydrogen  atom  at  1064  nm  for  various  laser  intensities 
(W/cm^),  compared  with  calculations  of  Krause,  et  al[l]  and  Potvliege  and  Shakeshaft  [2], 

Figure  2.  Same  as  Fig.  1 . 
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